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Abstract 


We  consider  an  unloaded  axially  aligned  collection  of  A  elastic  spheres  confined  within  rigid 
walls,  where  N  is  allowed  to  vary  between  several  and  several  hundred.  Alignments  with  both 
monosized  and  progressively  shrinking  spheres  (“tapered  chains”)  are  considered.  Being 
macroscopic,  the  spheres  are  assumed  to  incur  restitutive  losses  when  they  interact  and  we  adopt 
the  well-established  Walton-Braun  scheme  to  model  this  restitution.  These  spheres  are  known  to 
repel  upon  contact  via  an  intrinsically  non-linear  repulsive  force  law  upon  contact  due  to  Hertz 
(Hertz  law).  We  show  that  when  this  system  is  driven  at  an  edge  sphere  with  some  force  F  and 
frequency  Q,  this  driven  dissipative  system  can  respond  via  continuous  sequences  of 
overcompression  followed  by  dilation.  The  frequency  of  this  “breathing”  process  can  be  tuned 
across  a  wide  range  of  values  by  varying  the  geometric  and  mechanical  properties  of  the  elastic 
objects,  the  tapering,  and  the  system  size  and  thus  allows  for  resonance-like  behavior  of  these 
nonlinear  systems.  We  shall  discuss  how  this  system  offers  the  potential  to  develop  sensors  that 
may  be  able  to  operate  across  a  wide  range  of  frequencies  and  could  play  a  role  in  energy 
harvesting  from  a  wide  spectrum  of  sources  such  as  ocean  waves,  wind  and  possibly  geothermal 
sources.  The  project  has  also  led  to  two  PhDs  being  granted,  one  to  Dr  Robert  Paul  Simion  and 
the  other  to  Dr  Diankang  Sun.  Dr  Sun’s  work  was  closely  related  to  Dr  Simion’s  focus  but 
explored  mathematical  issues  related  to  grain  geometry.  In  addition,  Aharon  Festinger  did  an 
undergraduate  honors  thesis  on  trying  to  simulate  how  pulses  penetrate  into  loaded  3D  beds,  an 
extension  of  Dr  Simion’s  problem  to  3D  but  at  the  pre -breathing  level. 
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1 .  Introduction 


Granular  assemblies  provide  an  example  of  media  characterized  by  nonlinearity  and  discreteness 
[1-3].  Granular  alignments,  such  as  one  of  elastic  spheres  provide  one  of  the  simplest  and  richest 
dynamical  systems  that  exhibit  a  wide  range  of  unexpected  dynamical  behavior  [1,4].  Hence 
these  systems  have  been  the  object  of  intensive  theoretical,  simulational  and  experimental 
research  since  1983  [1,4]. 

Impulse  propagation  through  granular  chains  made  of  spheres  varying  in  size,  mass,  stiffness 
continues  to  be  a  heavily  investigated  subject  in  the  context  of  nonlinear  dynamics  [1,4]. 
Numerous  studies  have  shown  that  alignments  consisting  of  identical  spherical  grains,  which 
initially  barely  touch  one  another,  are  capable  of  sustaining  the  propagation  of  non-dispersive 
lumps  of  energy,  or  solitary  waves,  when  one  of  the  edge  grains  is  impacted  by  a  striker  [5-18]. 
If  the  mass  of  the  striker  has  parameters  (size,  density)  comparable  to  those  of  the  impacted  edge 
grain  a  single  solitary  wave  will  form  in  the  chain.  However,  if  the  mass  of  the  striker  is  much 
greater  than  the  mass  of  the  grain  then  multiple  solitary  waves  with  progressively  decaying 
amplitude  may  eventually  form  in  the  chain  out  of  an  initial  solitary  wave  train.  Sokolow,  Bittle 
and  Sen  [19]  found  numerically  that  in  a  granular  alignment  of  unloaded,  monodispersed  elastic 
spheres  of  mass  m  and  radius  R  a  solitary  wave  train  with  P  peaks  is  formed,  when  S  function 
perturbations  are  applied  to  P  grains.  This  work  has  since  been  experimentally  validated  by  Job 
et  al.  [20]  More  recent  studies  reveal  remarkable  shock  absorption  properties  of  small  horizontal 
granular  chains  made  of  progressively  decreasing  spherical  grains  [21]. 

A  granular  chain  in  its  simplest  form  consists  of  two  elastic  spheres  of  radii  Rt  and  Ri+l 
subjected  to  an  axial  compression  (without  exceeding  the  elastic  limits  of  the  constituent 
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materials).  The  repulsive  Hertz  potential  [22],  i.e.,  the  potential  energy  stored  at  the  contact 
between  the  two  elastically  overlapping  spheres -can  be  written  as: 


m,M) = 


RR 


J  /+1 


5  A,f+1 


R,  +  R 


s5/2 


i+l 


(1) 


where  Di  i+l  is  a  material  dependent  coefficient  and  5t  i+l  represents  the  overlap  between  the  two 
spheres: 

S,,M  =  Ri  +  -  (XM  ~xt)^°-  (2) 


When  both  the  grains  are  made  of  the  same  material,  the  factor  Z)  i+l  depends  on  the  Poisson’s 


^  j  _ (j  2 

ratio  cr,  and  the  bulk  Young’s  modulus,  Y,  and  D  =  —  ( — — — )  and  the  derivative  of  the  Hertz 
potential  with  respect  to  8 ,  the  force  felt  at  the  interface  between  the  two  elastic  beads,  is 


nsUi+ 0  = 


1  Ri^i+l  e  3/2 


D  v  R  +R 


03/  z  _  o  3/2 

^i,i+ 1  i,i+\ ) 


(3) 


when  the  spheres  are  in  contact  with  one  another.  If  the  spheres  lose  contact,  the  Hertz  repulsion 
is  no  longer  present  [22]  and  the  interaction  force  vanishes.  For  a  weakly  pre-compressed 
granular  chain  made  of  identical  spheres  in  which  a  solitary  wave  has  been  initiated  at  one  of  its 
extremities  via  a  short-lived  perturbation  -  an  impulse  function,  the  speed  of  sound  scales  with 
Fl/6  [4]. 

While  most  of  the  studies,  have  been  focused  on  understanding  the  impulse  propagation 
through  impacted  granular  chains  via  solitary  waves,  the  driven  granular  systems,  which  have 
not  yet  been  systematically  probed,  seem  to  present  interesting  nonlinear  features  leading  to  what 
we  call  nonlinear  granular  breathing  [23].  Granular  breathing  and  its  possible  consequences  in 
the  context  of  energy  harvesting  was  investigated  for  chains  constituted  of  identical 
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(monodispersed)  and  progressively  shrinking  (tapered)  spheres,  when  a  constant  force  F is 


applied  to  one  of  the  edge  grains  [23].  We  assume  that  F  is  directed  into  the  chain.  For  lossless 
monodispersed  and  tapered  chains  the  product  TFU(\  which  is  approximately  proportional  to  the 
length  of  the  chain  reads  : 


TF'16  oc  R 


4/3 

large 


l-(l-g)° 

l-(l-g)' 


1/2  rU/3 


pulD 


(4) 


with  a  «  1/3  (a  numerically  obtained  result)  where  T  is  the  breathing  period  -  the  time  elapsed 
during  one  energy  cycle,  F  is  the  magnitude  of  the  applied  force  that  acts  on  the  first  grain  at  all 
times  and  the  tapering  percentage  q  is  defined  as  q  =  1  -  Ri+l  /  Ri .  In  earlier  work,  Sen  et  al. 
[24]  have  suggested  that  by  coupling  with  appropriate  sensors,  these  nonlinear  features  could  be 
exploited  in  recovering/hamessing  some  of  the  energy  released  by  the  ocean  waves  as  they  come 
ashore. 

In  the  present  paper  we  analyze  the  dynamical  behavior  of  driven  granular  alignments  when  an 
external  force  varying  sinusoidally  with  time,  F  =  F0  sinQ? ,  is  applied  to  one  of  the  edge  grains 
(for  monodispersed  chains)  or  to  the  largest  grain  (for  tapered  chains).  We  present  details  about 
our  simulations  here  on  the  monodispersed  system  and  present  the  results  of  the  tapered  chain 
system  without  further  elaboration  on  the  tapered  chain  system  itself. 

We  look  for  the  “resonance-like”  phenomena  in  these  highly  nonlinear  systems  in  the 
presence  of  some  dissipation.  The  reason  why  the  sought  out  process  differs  from  resonance 
itself  is  as  follows.  Unlike  in  normal  resonance,  the  system’s  response  in  terms  of  how  energy 
varies  in  time  (recall  this  is  a  dissipative  system)  is  not  purely  oscillatory  but  rather  unique  in 
nature,  not  describable  by  any  simple  known  function  that  we  are  aware  of.  In  addition,  unlike  in 
linear  systems  that  are  characterized  by  harmonic  interactions  between  the  particles,  the  system’s 
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response  frequency  depends  on  the  chain  geometry,  the  material  medium  and  how  the  amplitude 


of  the  external  driving  force  could  influence  the  behavior  of  the  system. 

The  restitutive  losses  are  described  via  the  quantity  w  as  follows  according  to  the  Walton- 
Braun  scheme  [25] 


F 


unloading 


F, 


- 1  -  w. 


loading 


(5) 


where  Funloading  <  Floading  and  where  Funloading 


denotes  the  force  between  two  adjacent  beads  as 


they  depart  one  from  another  and  Floading  is  the  force  between  the  two  grains  as  they  press  one 
against  another  [25].  When  the  equality  holds  true  there  is  no  dissipation  in  the  system  while  if 
F unloading  <  Fhading  the  energy  is  primarily  lost  during  the  loading  phase.  For  a  certain  frequency  of 

excitation  at  or  near  one  of  the  natural  frequencies  of  the  system  damping  is  of  primary 
importance, 


2.  Newtonian  dynamics  simulations 

Our  systems  are  modeled  as  horizontal  chains  of  identical  or  progressively  shrinking  spheres 
with  their  centres  axially  aligned,  made  of  the  same  material  type  and  placed  between  two  rigid 
flat  boundaries.  At  the  instant  t  =  0,  when  the  spheres  barely  touch  one  another  (no  Hertz 
repulsive  interaction),  a  time-dependent  acceleration  is  imparted  to  one  of  the  edge  grains.  The 
nearest  neighbor  interaction  of  the  N  -  2  identical  spherical  grains  (without  the  edge  grains)  is 
described  by  the  following  coupled  nonlinear  equations: 


mixi=-(ai 


e3/  2 


A3/2' 

ai,i+l°i,i+l. 


i  =  2,...,N  —  1 


(6) 
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where  mi  is  the  mass  of  the  sphere,  xi  denotes  the  position  of  the  center  of  mass  of  grain  i,  and 
at  is  an  adjustable  parameter  which  depends  on  the  geometrical  and  material  properties  of  the 
grains.  As  initial  conditions  we  consider: 

xt  (t  =  0)  =  0,  xt  (t  —  0)  —  0,  i  -  1, . . . ,  N 

Fn  sin  (Q  t) 

xl(t  =  0)  = - ,  Xj(t  =  0)  =  0,  j  =  2,...,N  (7) 

mx 

Ri+ 1  =R^~q)  i  -  h--N 

where  t  e[0,?max]  denotes  the  time  variable,  tnvdx  sets  the  maximum  length  of  our  mns,  F0  is  the 

magnitude  of  the  applied  force  ,  Q  represents  the  driving  frequency,  m  the  mass  of  the  sphere 
and  the  dots  over  a  variable  means  its  time  derivatives.  R{  is  the  radius  of  the  largest  sphere  . 

The  equations  of  motion  are  numerically  solved  via  the  velocity  Verlet  algorithm  [26].  The 
position,  velocity  and  acceleration  of  each  grain  in  the  system  are  recorded  thus  the  kinetic, 
potential  and  total  energy  of  the  system  are  obtained.  We  have  used  different  time  steps  of 
integration  and  settled  for  A?  =  10“5  ps,  as  to  make  a  compromise  between  accuracy  and 
efficiency  and  our  runs  extended  over  104  ps.  Different  external  sinusoidal  forces  were  applied 
to  different  chains  with  different  geometries.  Metallic  grains  with  low  dissipation  coefficients 
have  been  the  object  of  numerous  experimental  studies  so  we  chose  T^AfiV  (p  =  4.42  mg/mm  , 
g-  0.34  and  E  =  llOGPa)  as  the  material  type  simulated  in  this  study  [27].  We  have  also 
simulated  and  obtained  similar  results  for  stainless  steel  and  SiC.  In  analyzing  small  systems, 
15  <  N  <  30  reasonable  values  have  been  assumed  for  the  restitution  coefficient  (  0  <  w  «  1)  and 
the  radii  of  the  grains  (5- 10mm)  such  that  the  results  of  this  study  could  be  compared  to  the 
results  obtained  from  a  real  experiment. 
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3.  Breathing  in  driven  chains  (same  material  medium) 

First  we  analyzed  the  dynamical  behavior  of  a  monodispersed  chain  (Rj  =  R,  mi  =  m ),  with  one 
of  the  outer  grains  subjected  to  an  external  time  dependent  force,  F  =  F0  sin  Q t . 

For  a  specific  chain  configuration  driven  with  a  force  of  magnitude  F0  we  tried  various 

driving  frequencies  Q  to  set  the  preliminary  stages  in  finding  the  nonlinear  resonance  design 
parameters  (number  of  grains,  particles  radii,  mass  densities,  and  material  type).  Phase  diagrams 
which  are  useful  tools  in  investigating  and  understanding  the  dynamics  of  nonlinear  systems 
were  the  starting  point  of  our  work.  In  a  Hertzian  system  (described  by  Eq.  (1))  a  significant  part 
of  the  input  energy  is  stored  as  potential  energy  at  the  grain-grain  and  grain-wall  contact 
surfaces.  Figure  1  shows  the  space-time  kinetic  energies  plots  for  a  system  driven  at  different 
forcing  frequencies.  The  chain  comprises  a  small  number  ( N  =  15)  of  identical  elastic  grains  and 
dissipative  effects  account  for  the  energy  lost.  In  each  plot  the  horizontal  axis  indicates  the 
particle  site,  the  vertical  axis  shows  the  time  elapsed  and  the  shading  gives  the  kinetic  energy  of 
each  particle.  The  darker  the  spot  the  more  kinetic  energy  the  corresponding  particle  has.  When 
the  energy  is  supplied  to  the  system  at  the  same  rate  as  the  one  as  at  it  is  being  lost  the  resonance 
phenomenon  occurs. 

One  can  recognize  that  not  all  of  the  plots  in  Figure  1  depict  an  ordered  pattern  of  the 
particles’  energies  -  a  behavior  expected  to  be  found  at  resonance.  The  only  energy  density  plot 
which  shows  a  repetitive  pattern  closer  to  the  expected  resonant  behavior  is  the  central  plot  in 
Figure  1.  It  can  clearly  be  seen  that  in  the  central  plot  the  particles  are  moving  in  a  certain 
fashion  as  opposed  to  the  other  plots  in  which  the  particles’  motion  is  somehow  hard  to  trace. 

Once  the  conditions  for  defining  the  resonance-like  behavior  in  our  nonlinear  systems  were 
set  we  plotted  in  Figure  2  the  temporal  behavior  of  the  system.  The  total  kinetic  and  potential 
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energies  resonate  at  different  frequencies,  which  are  results  of  the  fact  that  a  granular  chain  is  not 
a  conservative  system,  energy  is  continually  exchanged  with  the  driving  mechanism,  and  energy 
is  being  transferred  to  the  medium  and  lost  via  restitutional  losses.  The  total  energy  of  the  system 
rises,  reaches  a  maximum  and  then  falls  in  a  repetitive  pattern  (the  total  energy  peaks  have 
almost  the  same  magnitude).  When  dealing  with  macroscopic  discrete  systems  we  do  not  expect 
to  have  only  an  exact  value  of  the  resonance  frequency.  However,  if  the  driving  frequency  does 
not  closely  match  the  required  resonance  frequency  (is  either  well  above  or  below)  the  system 
does  not  experience  a  resonant  behavior,  as  shown  in  Figure  2. 

The  next  question  we  asked  was  how  the  driving  frequency  relates  to  N ,  the  number  of 
spheres  making  up  the  chain,  when  all  the  other  parameters  (F0  ,R,q,w,  p,D )  are  kept  constant. 

Our  numerical  results  show  that  when  only  the  number  of  the  grains  changes  the  driving 
frequency  changes  according  to  the  following  relation: 

£\NX  =D.2N2  and  q  oc  — 

N  (7) 


We  have  also  investigated  the  relation  between  the  driving  frequency  Q  and  F0 ,  the 
magnitude  of  the  driving  force  for  a  constant  configuration  (N  ,  R  ,  q  ,  w  ,  p  ,  D  )  .  Our 
simulational  results  show  (Figure  3)  that  Q  changes  with  F0  according  to  the  following  relation: 


F, 


1/6 


01 


F, 


1/6 


02 


Q, 


Q 


and 


F, 


1/6 


Q 


=  const  ,  FI  cc  F 


1  /  6 


(8) 


We  have  found  that  when  only  F0  and  N  vary,  the  following  relation  holds  true: 
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The  energy  transport  through  a  chain  of  equal  sized  spheres  depends  highly  on  the  radii  of  the 
grains.  If  only  the  radius  is  the  adjustable  parameter  then: 


Q, 

Q2 
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*1 


4/3 


and  Q  oc 
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R 4/3  . 


(10) 


By  using  all  the  above  results,  we  find  that  for  a  granular  alignment  of  identical  spheres  (same 
size  and  same  material  type)  subjected  to  a  time-dependent  force,  the  driving  resonance 
frequency  reads: 


F 

Q  oc  0 


1/6 


NR 


4/3 


(11) 


One  would  certainly  expect  that  Q  depends  on  the  characteristic  properties  of  the  material 
medium  ( D,p )  to  transfer  the  energy.  The  energy  transport  depends  on  both  on  the  inertial 


property  of  the  medium  (to  store  kinetic  energy)  and  the  elastic  property  of  the  medium  (to  store 
potential  energy).  If  the  spheres  are  made  of  the  same  material  type,  we  can  guess  that  the 
inertial  property  is  the  mass  density  p  .  For  the  elastic  property  we  can  assume  that  the  pre-factor 


D ,  which  contains  the  Young’s  modulus  and  the  Poisson’s  ratio,  would  be  a  suitable  candidate. 
Mathematically  this  could  be  expressed  in  he  following  form: 


F, 


1  /  6 


Q.NR 


4  /  3 


f(D,p,w) 


(12) 


where  the  function  f(D,p,w )  depends  on  the  material  medium.  Using  dimensional  analysis  we 
conclude  that: 


p  1/6 

°d4/3  =  const{D,p,w )  =  f(w,p)Dl/ 3 

1  INK 

and  continuing  our  dimensional  analysis  we  arrive  to  the  following  relation: 
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(13) 
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=  k(w)D  p 


where  k(w)  is  a  dimensionless  coefficient  that  depends  only  on  the  restitution  coefficient  w . 
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4.  Breathing  in  driven  chains  (different  material  medium) 

For  two  different  simulated  materials  (Ti6Al4V  and  stainless  steel)  with  the  same  force  amplitude 
(  F0 ) ,  and  the  same  N  having  the  same  radii  (R)  the  following  relation  applies: 


Z)“1/3 

J  Ti6Al4V 

We  have  carried  out  our  extensive  dynamical  studies  in  the  presence  of  restitutive  effects  with 
the  restitution  coefficient  (w  =  0.01)  defined  as  in  Eq.  (5).  From  our  simulations  it  could  be 
inferred  that  when  two  monodisperse  granular  chains  are  made  of  different  materials: 


77  i;6  -i/2 

QNR4/3 


77  1/(1  -!/2 

ro  P  Z)~1/3 

QNR 4/3 


=  k(w) 


where  k(w)  is  a  constant  for  a  given  w  .We  found  k  «  3.685  for  w  =  0.01.  To  clarify  how  a 


change  in  w  could  influence  the  resonance  behavior  more  studies  are  under  way.  At  present  our 
investigations  do  not  reveal  any  problems  regarding  the  validity  of  Eq.  (15)  for  values  of  w 
slightly  greater  than  0.01. 

For  a  monodispersed  chain  obeying  to  the  conditions  stated  above  our  studies  indicate  that 
the  driving  (resonance)  frequency  is  of  the  form: 


Q  = 


F 1/6  1 

_ j_o _ 1 

NR  413 pll2D 1/3  k(w) 


(16) 


We  mentioned  before  that  if  the  frequency  of  excitation  is  either  above  or  below  the  resonance 
frequency,  we  expect  that  the  system  exhibit  a  somewhat  chaotic  response  pattern.  For  a  linear 
system,  when  the  external  force  has  a  certain  forcing  frequency  of  oscillation,  the  steady-state 
response  will  have  the  same  frequency  of  oscillation.  However  a  nonlinear  system  could 
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presumably  exhibit  chaotic  behavior  at  frequencies  below  and  above  the  resonance-like 
frequency. 

If  the  magnitude  of  the  external  force  and  the  topology  (size,  material  type)  of  the  system  are 
not  altered,  we  find  that  more  than  one  driving  frequencies  can  induce  the  expected  resonance 
behavior  in  a  granular  chain. 

The  panels  in  Figures  4  and  5  show  that  only  at  or  near  certain  values  of  the  driving 
frequency  the  external  source  supplies  the  system  with  energy  at  a  rate  almost  equal  to  that  lost 
through  dissipative  effects.  The  shape  of  the  total  energy  peaks  change  with  Q  in  an  unexpected 
manner.  Note  that  these  results  are  not  numerical  artifacts  but  are  due  to  the  discreteness  and 
nonlinearity  of  these  systems  and  we  are  not  in  a  position  now  to  analytically  describe  the 
nonlinear  breathing  processes  described  in  this  work. 

We  find  these  unexpected  results  only  for  or  near  certain  driving  frequencies,  and  as  we 
lower  the  frequency  the  number  of  sub-peaks  in  each  peak  in  the  energy  versus  time  plot 
increases.  We  tried  to  find  a  relation  between  the  number  of  peaks  q  and  the  resonance 
frequency  Q  and  in  order  to  simplify  the  analysis  we  assumed  the  simple  form: 


Q  = 


p  i  i 

Kl  f(q,N)  ’ 


(17) 


where  (3  =  F01/6  / pV2Dv 3 


is  a  constant,  k(r/,w)  is  a  function  that  depends  only  on  the  number  of 


the  peaks  for  a  given  value  of  w-  0.01  .The  function  f(q,N) 


\-{\~q)Nn 
i  -  (i -q)in 


which  from 


our  previous  results  is  a  constant  for  a  given  configuration  (N,  q)  turns  to  N  in  the  limit  q  — »  0 
(  monodispersed  chain)  and  we  recover  Eq.  (16)  with  a  slightly  different  value  of  k . 

Our  approach  was  to  first  calculate  the  approximate  value  of  Q  for  a  monodispersed  chain 
by  using  Eq  (16),  then  to  tune  /  adjust  the  computed  value  to  the  resonance  frequencies 
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corresponding  to  different  tapered  chains  with  the  same  N  (between  15  and  20  )  but  different  q 


values  (between  4  and  9  %)  and  finally  compute  different  values  of  kip)  for  different  77  ’  s.  By 


averaging  the  values  of  k  over  the  number  of  our  trials  and  plotting  the  obtained  average  values 
of  k  versus  the  number  of  peaks  77  we  finally  arrived  at  the  following  linear  expression: 
k(q)  =  mq  +  b  ?  (18) 

with  the  integer  77  >  1,  the  slope  m  =  2.5705  and  the  intercept  b-  1.  Systematic  computational 


studies  reveal  that  the  resonance-like  process  in  ID  chains  with  one  of  the  edge  grains  subjected 
to  a  sinusoidally  varying  force,  could  be  attained  for  a  driving  frequency  of  the  following  form: 


F0l/6  1  1  1-(1  -q)1'3 

F  iarge  k  (t]  ,  w)  p^D"3  1  -(  1  -</)*' 3 


(19) 


In  order  to  test  the  accuracy  of  our  results  we  have  used  Eq.  (19)  for  different  chain 
geometries  (N,RlaTge,q)  and  different  force  magnitudes  and  found  that  the  relative  percentage 

error  for  Q  is  less  than  3%  for  small  systems  (15<tV<50)  and  less  than  6  %  for  large  systems 
(100  <  N  <  500).  We  suspect  that  this  disagreement  could  mostly  be  due  to  the  way  the 
restitution  coefficient  was  defined  in  our  computational  studies  and  different  approaches  are 
under  investigation. 

Equation  (19)  shows  that  for  a  chosen  configuration  of  a  granular  chain  we  can  induce  the 
resonance-like  phenomenon  by  only  modifying  the  amplitude  of  the  driving  force  with  the 
corresponding  frequency.  In  principle,  Eq  (16)  gives  us  a  different  kind  of  freedom,  one  in  which 
the  resonance-like  frequency  can  be  changed  by  adding  or  removing  one  or  more  of  the  grains 
from  the  system  Thus  these  discrete  systems,  which  are  intrinsically  nonlinear  and  hence  highly 
sensitive  to  external  perturbations  could  be  tuned  to  respond  in  a  great  many  ways.  This 
enormous  parameter  space  across  which  the  resonance-like  process  can  be  realized  suggests 
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that  these  systems  could  in  principle  serve  as  versatile  sensors.  Our  preliminary  analyses  suggest 
that  frequencies  between  0.1  Hz  and  104  Hz  can  potentially  be  realized.  Perhaps  an  even  broader 
range  is  possible.  To  make  the  transition  from  these  challenging  mechanical  systems  to  perhaps 
commercially  available  sensors  will  not  be  easy.  We  comment  below  on  how  the  work  presented 
here  could  potentially  be  useful  in  the  context  of  energy  harvesting  applications.  These 
comments  are  based  upon  our  work  in  progress,  which  will  be  detailed  in  future  journal 
publications. 

5.  Why  is  breathing  of  interest  in  energy  harvesting  applications? 

Imagine  an  array  of  cylinders,  each  containing  a  system  that  is  capable  of  breathing,  and  the 
array  to  be  buried  just  under  the  surface  of  a  beach  such  that  one  face  of  this  array  is  exposed  to 
the  ocean  waves  as  they  break  ashore.  It  may  be  possible  to  arrange  this  system  to  undergo 
sustainable  breathing.  If  there  are  piezoelectric  sensors  panels  at  the  surface  of  the  cylinder  array 
that  is  opposite  to  the  ones  exposed  to  the  waves,  then  the  mechanical  energy  of  the  waves  can 
not  only  be  efficiently  absorbed,  but  transported  into  voltage  drops  arising  due  to  the  response  of 
the  piezoelectric  sensor  panel.  Using  a  circuit  such  as  a  half-bridge  power  harvester  with  a 
leakage  resistance,  it  may  be  possible  to  feed  the  voltage  drop  onto  a  power  grid  [24].  Thus,  the 
energy  of  the  ocean  waves  crashing  ashore  across  large  distances  of  the  world’s  beaches  could  be 
harvestable  using  few  moving  parts.  If  we  note  that  the  power  in  the  waves  is  estimated  to  be  1- 
10TW  and  that  this  is  in  the  same  ballpark  as  the  world’s  power  needs,  harvesting  the  energy  of 
dissipating  waves  could  be  a  desirable  source  of  energy  [24].  However,  one  may  be  able  to  do 
better  -  because  the  complicated  mechanical  systems  may  actually  be  replaceable  by  fully 
electronic  systems  [24].  The  following  closing  section  elaborates  on  this  point. 
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Nonlinear  dynamics  has  begun  to  take  shape  as  a  research  field  in  relatively  recent  times  [4]. 
While  a  great  many  simple  systems  are  nonlinear,  isolating  the  effects  of  nonlinearities  to 
understand  and  exploit  its  consequences  is  another  matter.  Hence,  there  is  a  strong  need  to  have 
available  experimental  systems  with  a  control  on  dissipation  where  the  effects  of  nonlinearities 
can  be  systematically  examined.  Newcomb  and  coworkers  have  pioneered  the  development  of 
circuits  with  low  power  requirements  that  can  behave  as  nonlinear  systems  [28].  While  the 
development  of  such  circuits  is  challenging  and  time  consuming,  the  potential  benefits  can  be 
many,  such  as  the  development  of  chips  that  can  effectively  behave  as  complex  nonlinear 
systems.  Thus,  given  some  simple  input,  these  chips  could  be  expected  to  output  what  a 
nonlinear  system  might  yield  as  a  solution. 

We  are  currently  studying  the  breathing  problem  with  the  focus  to  realize  the  lowest  and 
highest  frequencies  and  the  range  in  between  using  Eq.  (19).  Our  objective  is  to  develop  circuits 
that  will  be  able  to  behave  effectively  as  these  nonlinear  breathing  systems.  If  successful,  we 
expect  these  circuits  to  find  many  applications  in  energy  harvesting  technologies  such  as  from 
wind,  ocean  waves,  geothermal  sources  and  more. 
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Figure  Captions 


Fig.  1  The  figure  shows  the  potential  (PE),  kinetic  (KE)  and  the  total  energy  (TE)  of  a  driven 
Ti6A14V  monodispersed  chain  versus  time.  Here  N=  15,  Riarge  =  5mm  and  Fo  =  100  N.  The 
inset  shows  the  PE  behavior  over  a  longer  time  span 

Fig.  2  The  figure  shows  the  behaviour  of  TF 1/6  versus  q  for  systems  of  two  different  sizes  (N=  15, 
20).  The  results  are  obtained  simulationally  for  F0  values  of  100  N,  200  N,  400  N,  and  800  N. 

The  plot  basically  shows  how  the  chain  length  shrinks  with  increasing  q  and  this  behavior 
agrees  with  the  form  shown  in  section  4. 1 

Fig.  3  Space  versus  time  plots  of  kinetic  energy  of  grains  (shown  in  grey  scale  with  darker 
implying  more  kinetic  energy)  for  monodispersed  (q  =  0%)  chain  of  N=15  Ti6A14V  beads  with 
R  =  5  mm.  Here  Fo  =  100  N,  w=  0.01  for  a  range  of  different  driving  frequency  (Q)  values 

Fig.  4  Total,  kinetic  and  potential  energy  versus  time  behaviour  in  monodispersed  chains  of 
N=15  Ti6A14V  beads  at  driving  frequencies  Q=  2570  rads/sec  which  is  below  (upper  panel)  and 
at  Q=  3420  rads/s  which  is  above  (lower  panel)  the  non-linear  resonance  frequency  of  Q=  2990 
rads/s 

Fig.  5  Non-linear  resonant  behaviour  of  monodispersed  tungsten  bead  chain  with  N=  15, 

R  =  5mm,  w=  0.01,  F  =  100N  and  Q=  2171  rads/s  and  q=l 

Fig.  6  Two  distinct  non-linear  resonant  behavior  patterns  corresponding  to  r\  =  2,  3  (upper,  lower 
panels)  for  beads  of  Ti6A14V  with  q  =  5%,  N=  15,  F  =  100  N,  R  =  5mm,  and  w=  0.01.  For  r\=2, 
Q=  1709  rads/s  (upper  panel)  and  for  q=3,  Q=  1200  rads/s  (lower  panel) 
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Tabic  1 

Definitions  of  mathematical  symbols  used 

Symbol 

Definition 

^avrragK 

D 

Fa 

P 

1  unloading 

*< 

T 

w 

Y 

Zi 

Average  energy  transport  speed  in  the  chain 

Amplitude  of  applied  force  on  the  first  grain 
Force  on  the  grains  during  loading 

Force  on  the  grains  during  unloading 

Taperi  ng  coefficient 

Radius  of  the  Ith  grain 

Round-trip  energy  transport  time  in  the  chain 
Restitution  coefficient 

Young's  modulus 

Position  of  the  grain  center  in  the  unperturbed 
state  of  the  chain 

&U  *  1 

n 

p 

IT 

n 

Overlap  between  adjacent  grains  i  and  i  +  l 
Number  of  peaks  in  each  breathing  cycle 
Density  of  elastic  grain 

Poisson  ratio 

Driving  and  non-linear  breathing  frequency 
at  non-linear  resonance 

23 


0.035 

0.03 

0.025 

0.02 

S 

ui 

0.015 

0.01 

0.005 

0 

Fig.  1 


0  100  200  300  400  500  600  700 

t(j»s) 


24 


1.« 


Fig.  2 


25 


t(|js)  t(Ms)  t(ps)  t(Ms)  t(MS)  t(Ms)  t(Ms) 


10000 

9000 

10000 

9000 

- 

8000 

8000 

L.  i  - 

7000 

f  - 

* 

7000 

6000 

/  •  *“  ^ 

I  ■ 

*  * 

6000 

_s 

5000 

•  .• 

L- 

i  x' 

5mn 

' 

4000 

r 

► 

l 

4000 

A* 

3000 

3000 

2000 

2000 

- 

1000 

1000 

- 

0 

i  i 

0 

— 1_ 1  A 

51015  51015 

N  N 


10000 

- T - 1 - 1 

10000 

9000 

9000 

- 

8000 

8000 

-  - 

7000 

ap r»  r 

7000 

6000 

6000 

5000 

5000 

" 

4000 

4000 

- 

3000 

2000 

3000 

2000 

- 

1000 

1000 

ol 

i  i 

0 

K  J  i 

51015  51015 

N  N 


10000 

9000 

8000 

7000 

6000 

5000 

4000 

3000 

2000 

1000 

0 


10000 

- 

sr-  • 

^  - 

9000 

> 

-  ’  r  r 

•;V 

8000 

Jj*  - 

7000 

9 

IS-  - 

6000 

- 

5000 

\  -  . 

7*  -J 

>< 

>« 

4000 

S 

a- 

Sire  "" 

3000 

pT3 

2000 

. 

1000 

0 

i i  j 

51015  51015 

N  N 


10000 
9000 
8000 
7000 
6000 
5000 
4000 
3000 
2000 
1000 
0 

51015 

N 


Fig.  3 


26 


0.015 


0.01 


UJ 


0.005 


1000  2000  3000  4000  5000  6000  7000  8000  9000 

t(fiS) 


0.025 


002 


0.015 


UJ 


0.01 


0.005 


1000  2000  3000  4000 

t(HS) 


5000  6000  7000  8000 


Fig.  4 


27 


Fig.  5 


28 


29 


